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I. Introduction 

This report describes the progress made in the first year of a two-year 
project which deals with modeling the spatial distributions of observed come- 
tary radicals. 

The primary goals of the research program are the inclusion and evalua- 
tion of the following processes which are usually neglected in modeling 
cometary atmospheres: (1) dependence on heliocentric distance of initial 

comet expansion velocities and nucleus vaporization rates, (2) dependence on 
heliocentric velocity of photochemical lifetimes and photon emission rates 
where applicable, (3) the correct isotropic ejection of daughter radicals upon 
photodissociation of parent molecules, (4) solar radiation pressure, and (5) 
collisional effects on radical kinematics. The principal results of the 
modeling effort would be to generate the information necessary to reduce 
observed radical column abundances to radical production rates for CN, C2, 

Cg, and OH. Secondary goals which would naturally follow would be a better 
quantitative understanding of the physical processes which are important in 
cometary atmospheres. In all instances model development is constrained by 
the observed spatial distributions of these radicals already published in the 
literature, as well as by the best currently available photochemical data. 

In addition, the same models developed here could also be easily applied to 
analyze the new observed spatial distributions of cometary radicals. 

A fairly large amount of the filter photometric data can be expected to 
be acquired during the 1985-1986 apparition of Comet P/Halley from the members 
of the Photometry/Polarimetry network organized under the International Halley 
Watch. The large amount of high quality data requires much more sophisticated 
models for analysis than are currently available. 

The type of model developed is a many particle-trajectory Monte Carlo 
model. The facility and economy of using this type of model for treating the 
complex geometry and dependencies mentioned above has been discussed and 
demonstrated in the published literature (Combi and Delserame 1980a and 1980b, 
Combi 1980). The traditional approach initiated by Eddington's (1910) foun- 
tain model and continued in such later works by Haser (1957, 1966), Wallace 
and Miller (1958), Keller and Meier (1976) and Festou (1981) becomes computa- 
tionally unmanageable for multi-dimensional time-dependent models. This is 
also true for the one-dimensional fluid models which have begun to model the 
physics correctly (Marconi and Mendis 1982, Huebner and Keady 1982). 
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During the first year of the project all of the new processes have been 
incorporated into the Monte Carlo particle-trajectory model (MCPTM). Prelimi- 
nary analyses as to their effects on the spatial distributions of observed 
cometary radicals have been begun. These first exploratory model runs will be 
discussed in detail in the following sections of this report. 

Also, in the interest of obtaining new data regarding the spatial distri- 
butions of neutral cometary species, a collaborative effort with Dr. Uwe Fink 
of the University of Arizona has been established. 
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II . Observations of the Spatial Distributions of Cometary Radicals 

An important part of this project is to assess the observed spatial dis- 
tributions of cometary radicals in order to constrain the parameters of the 
MCPTM. Most observations of this type have been analyzed in terms of Haser's 
(1957) model. Although Haser model scale lengths cannot be used to give 
direct physical quantities (lifetimes), they can be useful tools to 
characterize the apparent size of the source and decay regions of cometary 
radicals. 

Newburn and Spinrad (1984) have presented Haser scale lengths for the 
parents of C 2 , C 3 and CN in several periodic comets. These were obtained by 
two-point filter photometry: one point centered on the photometric nucleus 

and another displaced some known distance from the nucleus (typically 1 to 
4 x 10^ km). Since Haser's model is basically a three parameter model (the 
production rate and two scale lengths), Newburn and Spinrad had to assume the 
values of A'Hearn (1982) for the radical decay length. From two points they 
could then calculate both a parent scale length and the production rate. 
Although errors could arise from both this assumption and that of spherical 
symmetry (i.e., no radiation pressure distortion), a large number of obser- 
vations may average out some of the random discrepancies. 

Figure 1 shows a plot of the CN parent scale length versus heliocentric 
distance for the combined data sets of Newburn and Spinrad (1984), Combi and 
Delsemme (1980) and Delsemrae and Combi (1983) as well as a newer value deter- 
mined by Johnson, Fink and Larson (1984) from their CCD long-slit spectra of 
Comet P/Tuttle. There is no gross systematic difference between the two data 
sets. The scale lengths of Newburn and Spinrad were determined from the two- 
point photometry of smaller short period comets at medium to large heliocen- 
tric distances, whereas our data and that of Johnson, Fink and Larson were 
determined from the average of the sunward and antisunward brightness profiles 
of larger new or very long period comets at small to medium heliocentric 
distances. There are large uncertainties involved in determining one scale 
length from only two points. Newburn and Spinrad admit a factor of 2 change 
in the parent scale length results from a 25% change in one of the two 
measured fluxes. On top of this is their need to assume A'Hearn's adopted 
daughter scale length law for CN, where we (Combi and Delsemme 1980b) had 
found in fact no heliocentric distance-dependent law for CN but a distribution 
of values varying rather randomly over an order of magnitude. Therefore, the 
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Figure 1. CN Parent Scale Length vs. Heliocentric Distance. The points 
plotted for Comets P/Encke , P/Tuttle, P/Stephan-Oterma and 
Bradfield were computed from two-point filter photometry by 
Newburn and Spinrad (1984). Those for Comets Bennett, Kohoutek 
and West were computed from entire brightness profiles (Combi and 
Delsemme 1980b, Delsemme and Combi 1983). The solid line is the 
best fit power law to the combined data sets. 
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facts that both sets of data indicate a similar dependence and the mean 
values agree to within a factor of 2, are quite encouraging though possibly 
just fortuitous. 

Taking all the data together, a power law can be fit which yields the 

straight line in Figure 1 and which has the form Ap = 1.6 x 10^ r^ 1,44 for \> 

in km and r^ in AU. The slope for this fit has an uncertainty of about ±0.3 

12 IS 

generally excluding both the typically adopted r^ or r H laws. An r H law, 

on the other hand, would be expected for a photodissociation lifetime 

2 “0 5 

(x « r^ ) and a variable parent velocity (v <* r^ ). 

Figure 2 shows the variation of the C 2 parent scale length with helio- 
centric distance for the combined data sets. Unlike the case for CN, though, 
there does seem to be a systematic difference. The scale lengths of Newburn 
and Spinrad tend to be both larger and exhibit a flatter slope than do those 
compiled by Combi and Delsemme (1985). Taken together, a power law in r^ can 

be fitted to all the data which have the form A_ = 2.2 x 10^ r H ^*^. The full 

P H 

profile data taken alone has the form \p = 1.6 x 10 4 r^* , whereas the data 
of Newburn and Spinrad yield A p = 3.8 x 10^ r^'^. 

At this point, there could be two reasons for such a difference. First, 
the C 2 spatial distribution in the smaller short period comets may be substan- 
tially different from that in the more active new and very long period comets. 
If C 2 were produced primarily by gas phase chemical reactions or an icy grain 
source as suggested by A'Hearn and Cowan (1980), there could certainly be 
differences between the two populations of comets. Second, there could be 
differences between the two model fitting procedures. The C 2 radical experi- 
ences nearly twice the acceleration and is likely moving at a lower speed than 
is CN (Combi and Delsemme 1985). Since the two-point photometry method does 
not average out the asymmetry, systematic differences could result. In addi- 
tion, the general problem discussed above where a modest uncertainty in one of 
the measured fluxes causes a large uncertainty in the determined scale length 

still applies. In any event, the C 2 results for complete brightness profiles 

0 

seem to present a reasonable case for an r^ law for the C 2 parent scale 
length. 

Finally, for Cg the only complete data are those by Newburn and Spinrad, 

o 

which also seem to be consistent with something like and r^ law for the 

o 2 

parent Haser scale length given by 2.5 x 10 r^ km. 
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Figure 2. C 2 Parent Scale Length vs. Heliocentric Distance. The points 
plotted for Comets P/Encke, P/Tuttle, P/Stephan-Oterma and 
P/Swif t-Gehrels were computed from two-point filter photometry by 
Newburn and Spinrad (1984). Those for Comets Burnham, Tago-Sato- 
Kosaka, Bennett and Kohoutek were computed from entire brightness 
profiles (Combi and Delsemme 1985). The solid line is the best 
fit power law to the combined data sets. 
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III. Photochemical Data 


Another important aspect of modelling the spatial distributions of come- 
tary radicals is the photochemical data relevant to the photodestruction 
(dissociation and/or ionization) of radicals and their prospective parents, as 
well as the fluorescence emission rate per indicent solar radiation flux (or 
g-factor). 


CN. 

Because of the existence of many solar absorption lines in the region of 
~ 3900A°, the g-factor for the CN(O-O) is not only a function of heliocentric 
distance (~ 1/r^ ) but is also a rather strong function of heliocentric velo- 
city (the Swings effect). This calculation has been performed by Tatum and 
Gillespie (1977) and Mumma et al. (1978). Both calculations demostrate the 
same relative heliocentric velocity dependence but differ by an overall factor 
of ~ 20% due to the respective use of disk center and disk average solar 
spectra (see A' Hearn 1982 for a plot of the two calculations). Figure 3 shows 
a plot of the calculation of Tatum and Gillespie adopted here as the only 
formally published data. The g-factor is not only important for calculating 
absolute column densities but also because it directly contributes to the 
radiation pressure acceleration which visibly distorts the CN coma (Combi 
1980). 

Our initial suggestion for HCN as the parent of CN (Combi and Delsemme 
1980b, Combi 1980) has remained a quite viable explanation of the observed 
brightness distribution. Cochran (1982) has reached a similar conclusion from 
observed brightness profiles of CN in Comet P/Stef an-Oterma. Huebner (1983) 
has recently revised the photochemical lifetime of HCN downward only slightly 
from 8.3 x 10^s (Huebner and Carpenter 1979) to 7.7 x 10^s, with the major 
contribution still from solar Lyman a. As before, an excess energy of 4.3eV 
divided among the fragments yields a random velocity for CN radicals of 
1.1 km/s with respect to the radially outflowing HCN molecules. 

c 2 . 

The g-factors for the various bands of the Swan system of C 2 are rela- 
tively well known and there seems to be no significant heliocentric velocity 
dependence (A' Hearn 1975, A' Hearn 1982, Danks and Lambert 1983). The case for 
identification of the production mechanism for C 2 (i.e., possible parent 
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HELIOCENTRIC VELOCITY (KM/S) 

Figure 3. Fluorescence Excitation Rate for CN(O-O) at 3883 A. Plotted is the excitation rate 
(g-f actor for the (0-0) band of CN at 3883 A as a function of heliocentric velocity 
as computed by Tatum and Gillespie (1977). 



molecules) is, however, not so advanced. A recent paper by Cochran (1985) 
has advanced the theory that the shape of three C 2 brightness profiles in 
three different comets can be explained by the photodissociation chain 
C 2 H 2 ♦ C 2 H -*■ C 2 provided that the unknown lifetime of C 2 H is about 2000 
seconds. This conclusion was based on a one-dimensional single-fluid chemical 
model (Cochran 1982) which included the estimation of the C 2 H 2 photodissoci- 
ation lifetime of 3.2 x 10^ seconds originally calculated by Huebner and 
Carpenter (1979). The conclusion then really implies that effective parent 
lifetime for production of the C 2 radical is of order 3.4 x 10^s and that the 
lifetime of C 2 H must be negligible ( 2000 s) as compared with that for C 2 H 2 
itself (~ 3.2 x 10^s). Very recently, however, Huebner (1985) has revised the 
early study of acetylene photochemistry and finds the following branches and 
rates are dominant: 



-1 

-5 -1 

+ hv + C 2 H + H 

T = 1 X 

10 s 


“I 

-6 -1 

* C 2 + H 2 

t = 2.7 

x 10 s 


Therefore, the new lifetime for acetylene is 7.9 x 10^s, which is more than a 
factor of two larger than that implied by the observed profile. This would 
seem to eliminate C 2 H 2 as the only or even primary parent for C 2 . 



Very little is known photocheraically about C^ or prospective parents. 

The g-f actor for the -*■ X 1 Z + band at 4040A 0 has recently been revised 

u g 

upward by a factor of 40 (A'Hearn 1982). Thus, instead of being as nearly 

abundant as C 2 and CN as was believed, production rates are now reduced to 

over an order of magnitude less. One implication of this is that C 3 can no 

longer contribute significantly as a source of C 2 by photodissociation. 

Another is that the radiation pressure acceleration of C 3 radicals by solar 

_o 

light can now be calculated to be .56 cm s , which is larger than both C 2 and 

—2 

CN, .42 and .34 cm s respectively. 

The simplest prospective parent for C 3 for whose photochemical lifetime 
has been calculated is CH 3 C 2 H (propyne). Stief (1972) has suggested the 
following reaction chain in which propyne is photodissociated leaving C 3 H 2 in 
an excited state that spontaneously decays to form C-^: 
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CH 2 C 2 H + hv C^H 2 + H 2 
C3H2 •+■ C3 + H2 

Stief found a rough lifetime for this reaction of 5000s, which is actually 
reasonably consistent with the Haser scale length of 2.5 x 10 km found by 
Newburn and Spinrad (1984, see Section II) for an expected parent velocity 
of 58 km s 

OH. 

Schleicher and A'Hearn (1982, 1984) have calculated the dependencies of 
both the OH photodissociation lifetime and the OH solar fluorescence rate on 
heliocentric velocity. At a distance of 1 AU from the sun, the OH lifetime 
can vary from £2 x 10^ s to ~4 x 10^ s depending on the heliocentric velo- 
city. The (0-0) band fluorescence rate varies by more than a factor of 6 for 
different heliocentric velocities. Figures 4 and 5 show these dependencies 
which have been incorporated in to the cometary radical particle-trajectory 

model. It is important to note that these calculations are made only for 

2 2 + 

transitions between the ground X n and the A £ electronic states. Fluor- 
escence occurs for transitions up to the v' = 0 and 1 vibrational levels, and 
predissociation occures for v'> 2. Furthermore, Singh, van Dishoeck and 
Dalgarno (1983) point out possible transitions to other electronic states with 
an important role for Lyman = absorption. Van Dishoeck and Dalgarno (1983) 
suggest a combined rate of 3.8 x 10 - ^s - ^ for these other states. , Simply 
adding the rates from the two calculations together would yield a total photo- 
dissociation lifetime in the range of only 1 to 1.6 x lO^s, also decreasing 
the magnitude of the heliocentric velocity dependence. However, Schleicher 
(1983) points out that Festou (1981b) using the vectorial model and the 
correct kinematics for ^0 photodissociation finds lifetimes for OH in the 
range of 2-3 x 10^s, apparently just in the range predicted by the Schleicher 
and A'Hearn calculation alone. 

IV. The New Monte Carlo Particle-Trajectory Model 

•The original Monte Carlo particle-trajectory model (MCPTM) as developed 
by Combi and Delsemme (1980a) has been generalized during the first year of 
this project and now includes: 
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(1) solar radiation pressure 

(2) isotropic ejection of daughter radicals owing to excess photolysis 
energy, 

(3) the heliocentric distance dependent parent velocity (Delsemme 1982), 

(4) time (heliocentric distance) dependent source rate and lifetimes, 

(5) heliocentric velocity dependent lifetimes and excitation rates, and, 

(6) multiple elastic scattering of cometary radicals by the outflowing 
g as * 

As each new process was incorporated into the MCPTM, exploratory model runs 
were made in order to test the effect of the process on the observed spatial 
distribution. I will now present separate discussions for each new physical 
process giving a brief description of the modeling procedure and some pre- 
liminary model results. 

A. Radiation Pressure and Variable Parent Outflow Velocity 
Although the original model included radiation pressure (also see Combi 
1980), a constant velocity for the outflowing parent molecules had been 
assumed. Furthermore, no calculations had ever been carried out to invert 
nucleus-centered filtered photometric observations with a model which included 
both the isotropic ejection of daughter radicals (also as in the vectorial 
model) and radiation pressure. As mentioned in Section II of this report, 
Haser parent scale lengths determined from observed CN profiles indicate some- 
thing like an r^*^ law. This would be expected for the photodissociation of 
parent molecules combined with the outflow velocity law determined semi- 
empirically by Delsemme (1982) to be v = .58 r^ ' km s . The production 
rates for CN as determined from the photometry of Comet West (1976VI) (A'Hearn 
et al. 1977, A'Hearn and Cowan 1980) have been recalculated for various model 
process combinations in order to test their effects. 

A'Hearn et al. (1977) had originally used an law for a Haser model 

parent scale length for CN and found that the production rate for CN varied 

_n Q 

with heliocentric distance as r^ . Combi and Delsemme (1980b) then 

measured CN parent scale lengths in Comets Bennett and West and found a de- 
1 8~h 1 

pendence of r^ * which they concluded was reasonably consistent with an 

o 2 

r H law. Both they and A'Hearn and Cowan (1980) then adopted the r^ law and 
found a CN production rate varying approcimately as r H _n where n “ 1.6±0.3. 

For this study, the CN production rates have been calculated using three 
different model descriptions and compared with the Haser model parameters 
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A' Hearn and Cowan have adopted. The results of these four cases are shown in 
Table 1. The first case is a Haser model, using the Average Random Walk Model 
(Combi and Delsemme 1980a) corrections to the scale lengths. Here the rjj 
variation in parent velocity and the isotropic ejection of the daughter 
radicals is taken into account in an approximate way. The second case is a 
Monte Carlo Particle Trajectory Model (MCPTM) which explicitly calculates the 
effects included in the first case. And the third case is a MCPTM which also 
includes radiation pressure (Combi 1980). In all cases the heliocentric 
velocity dependence of the CN(O-O) band excitation (Tatum and Gillespie 1977) 
is included. 

The largest effect is that of the variation of the parent velocity which 
steepens the slope of the resulting production law and is present in all three 
cases. The explicit inclusion of isotropic ejection (case 2) and radiation 
pressure (case 3) each flatten the slope somewhat but are nonetheless impor- 
tant. A comparison of the original Haser model calculations and the full 
MCPTM (Case 3) for the CN production in Comet West is shown in Figure 6. 

B. Multiple Elastic Scattering of Cometary Radicals 

Models of neutral cometary comae have generally dealt with the two 
extreme cases: one where the densities are high enough so that all neutral 

species are thermalized and are characterized by a bulk radial flow away from 
the nucleus (e.g., one-dimensional fluid models by Giguere and Huebner 1978, 
Marconi and Mendis 1982, Cochran 1982), and the other where the densities are 
so low that essentially no collisions occur (e.g., the vectorial model by 
Feston 1981a). 

A characteristic distance from the nucleus separating these two regions 
has been defined where the mean free path for a collision is equal to the 
distance from the nucleus. Whipple and Huebner (1976), for example, express 
this collision radius as 


r 

c 


= _ 2 £ 
4uv 


where Q = molecular production rate 

— IS 2 

a = cross-section for elastic scattering ~10 cm 
v = molecular outflow velocity. 
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Table 1 


CN Production in Comet West 




Model 


Ql ( s 

l)( a) 

(b) 

n 


A’ 

Hearn and Cowan Haser Model 

1 . 1 X 

10 27 

-1.61 

1) 


Average Random Walk Model 


1 .4 x 

10 27 

-2.00 

2) 

Monte 

Carlo Particle Trajectory 
(no radiation pressure) 

Model 

1.2 x 

10 27 

-1.92 

3) 

Monte 

Carlo Particle Trajectory 
(with radiation pressure) 

Model 

1. 3 x 

10 27 

-1.87 


Production Rate at 1 AU 
Exponent in power law Q = Qj r H ' 
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Figure 6. CN Production Rate vs. Heliocentric Distance. The CN Production rates 
for Comet West are shown as computed both with the Haser model scale 
lengths .of A'Hearn and Cowan (1980) 1 0 ] and with the Monte Carlo 
particle trajectory model which includes radiation pressure and a vari- 
able parent velocity [.+] . Intermediate cases are discussed in the text. 



The value of r c varies from ~10^ km for a bright comet like Kohoutek near 
perihelion, down to ~10^ km for a small short period comet like P/Encke at 
1 AU. 

In reality, though, there is no clear separation between the two zones. 

In order for a typical radical (CN, C 2 , OH) which has been emitted from its 
parent with some exothermic velocity to be truly therraalized, it needs on the 
order of 10 collisions with outflowing water molecules. On the other hand, 
many radicals produced outside of r £ will be subjected to at least 1 to 5 
collisions, thus questioning the true applicability of a free flow model. 

What we truly have, then, is a large regime on the order of 0.2 r c to 5 r c 
where partial thermalization will occur. 

The true random walk nature of the elastic scattering problem is modeled 
in a straightforward manner in the Monte Carlo particle-trajectory model 
(MCPTM). Kitamura et al. (1983) have included collisions in a Monte Carlo 
model essentially based on our original model (Combi and Delsemme 1980a). 

Since they basically followed the correct procedure, we have adopted a similar 
algorithm. There is one problem with their procedure which involves the ex- 
pression they assume for the density of outflowing molecules with which the 
modeled neutral radicals are colliding. They assume a Haser-like description 
for the density which is 


where Q, v are as above 

r = distance from the nucleus 
TH 2 O = photochemical lifetime for ^0. 


Strictly speaking, this is the correct density distribution for ^0 molecules, 
but as the H 2 O molecules dissociate, they produce OH and H, which still 
provide targets with which radicals collide. Since the effective cross- 
sectional area for the fragments is larger than the parent H 2 O but their 
outflow velocity is, on the average, also larger, a better approximation would 
be to simply drop the exponential term in the density expression. This yields 


H 2 0 


(r) 


Q 

Atrvr ^ . 
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For the simple case of no radiation pressure, this yields a closed form 
algebraic expression for the collision length of an arbitrarily moving neutral 
species, in place of the integral found by Kitamura et al. For the radiation 
pressure case the parabolic trajectory is divided into finite elements each of 
which is treated as a straight-line case. 

The random walk algorithm begins with the production of a neutral radical 
by its parent with a given velocity vector. From this location and velocity 
vector, a collision length is given by assumed values for the collision cross 
section, the molecular (l^O) production rate, and a random number. The 
neutral is then displaced this distance, scattered elastically off a moving 
water molecule and continued on for another collision. The trajectory 
proceeds this way until the observation time snap-shot. 

Preliminary model runs have been performed for the case of CN produced by 
the photodissociation of HCN (Combi and Delserame 1980b). We have chosen para- 
meters appropriate for the photochemical lifetimes at 1 AU and have neglected 
radiation pressure, which is of less importance at 1 AU, for the purpose of 
simplifying the preliminary analysis. The physical parameters of the coma 
models are: 


r H = 1 AU 
x(HCN) = 9 x 10 4 s 
t(CN) = 2.75 x 10 5 s 
v(parent) = 0.58 km/s 
v e (CN) = 1.02 km/s 

Scattering = 1 x 10 " 15 cm 


-2 


Models were run for cases of general molecular production rates from 
1 x 10 s to 1 x 10 s . The collision statistics from these model runs 
summarized in Table 2. As can be seen, for larger gas production rates, only 
half of the CN radicals produced are in a true free-flow regime. Furthermore, 
the effect collisions have on the determination of radical production rates is 
also evident in these model runs. For the case of the largest gas production 
rate, if one neglects collisions in the model, an underestimate of 11% to 15% 
in the determination of radical production rate would be made, depending upon 
the size of the photometer aperture used. The systematic errors, when re- 
ducing data for even smaller heliocentric distances (r H ~ 0.5 AU), should be 
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Table 2 


Collision Statistics for the CN Model at 1 AU 


Fraction of Radicals Undergoing n Collisions 


Qh 2° 

0 

1-2 

3-5 

6-10 

>10 

1 x 10 28 

.723 

.272 

.0042 

.0005 

.0003 

5 x 10 28 

.705 

.283 

.0094 

.0023 

.0005 

1 x 10 29 

.681 

.296 

.0160 

.0055 

.0009 

2 x 10 29 

.653 

.311 

.0247 

.0101 

.0016 

5 x 10 29 

.576 

.343 

.0503 

.0252 

.0055 

1 x 10 30 

.503 

.355 

.0852 

.0459 

.0147 


19 


quite important in bright comets such as P/Halley, since at these distances 
the collisional effects are compounded by larger gas production rates and 
smaller radical parent scale lengths. 

C. Time Dependence 

Time dependence in comet model parameters such as production rate, 
photochemical lifetimes and velocities enters naturally through the changing 
of the comet's heliocentric distance and velocity along its orbit with time. 
Models by Keller and Meier (1976) and Cucchiaro and Malaise (1982), for 
example, have treated time-dependent production rates, with the former also 
treating the time-dependent hydrogen lifetime and radiation pressure accelera- 
tion. As part of this project, however, time dependence has been included in 
all aspects of the extended source, kinematics, and decay of observed cometary 
radicals. The ability to handle such a general calculation illustrates the 
power of particle-trajectory models using Monte Carlo techniques. 

The original calculation, described in detail by Combi and Delsemme 
(1980a), and since then adopted by others (see Bockelee-Moran and Gerard 1984, 
Kitamura, Ashihara and Yamamoto 1985, and Schloerb and Gerard 1985), corres- 
pond to steady-state production rates and constant photochemical lifetimes for 
parent (and grandparent) and daughter species. Typically, N (~ 10 5 ) parent 
molecules are emitted at random times between the current or observation time 
and some back-up time tg which is long enough to build up the whole observed 
cloud. In this case, the model production rate is given simply by N/tg, which 
may be simply scaled to the desired production rate or fitted to an observa- 
tion to yield the production rate. Given typical random numbers (R^) on 
the interval [0,1) the back-up time for emission of a parent is given by 

C I = R i t B* 

In the steady-state MCPTM the manner in which the constant lifetimes for 
both parent and daughter decay has been modified from the original model in 
the interest of improving the statistical coverage. In the old procedure a 
decay time was computed using the principal of Monte Carlo, which yields the 
relation t Q = -t In (1-R^) where t is the exponential decay lifetime. This 
was computed twice, once for the parent and once for the daughter; a daughter 
trajectory was only computed only if the parent decay occurred within a time 
t I anc * the daughter decay did not. A better method is to use the concept of a 
forced dissociation for the parent and a simple weighting function for the 
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daughter decay (see forced first scattering discussion by Cashwell and Everett 
1959). A parent dissociation within the time t^ is given by 

t = -r In • (l - R. [l - e” tl/T P]) 

D p x L 

and the daughter is assigned a weight (* 1) given by 

w - p - r c i /T p) 


where Xp and x^ are the parent and daughter decay lifetimes. The first term 
(in parentheses) results simply from the forced parent dissociation. In the 
steady-state MCPTM the trajectories are calculated and collected in space and 
column density bins as before. 

In the time-dependent model, the variable production rate is treated 
quite simply by weighting a parent molecule according to the production rate 
at the emission time ^observation - relative to the production rate at 

the observation time. Thus a water vaporization curve or a simple power law 
in heliocentric distance can be assumed; for that matter, transient activity 
or outbursts can also be accounted for by simple weighting. 

The variable lifetime for the parent and daughter decay are somewhat more 
involved. The generally time-dependent lifetimes x = x(t) can be approximated 
as a set of finite element steps 


r [r^(T^), X2(T£ ) , x^T^), r n (T^)] 


where (T^, T 2 , T 3 ,...T n ) are sub-intervals in time such that 


T. + T-, + To + ... + T = t T . 
1 i. i n i 


The partial weight associated with the forced dissociation is then given by 


W Q = 1 - exp[-(T i /x i + T 2 /x 2 + T 3 /x 3 + ... + T n /x n )] 


The probability for a dissociation somewhere on this interval 0 < t^ < t^ is 


P(t D ) = 1 - e p 
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where 



T 

m- 


1 


T 

m- 


1 


+ t 


D 


(T. + T„ + ... + T ) 
i L m 


T 

m 


and 


T, + T„ + ... + T < t_ < T, + T„ + ... + T +...+T. 
12 m D 1 2 m n 


The quantity, p, is just the number of mean free times (like mean free paths) 
contained in time tjj. 

Just as in the steady-state case, the random number enters as the 
fraction of parents remaining after the time t^; therefore, 


R ± - (1 - e" p )/W D 


and solving for p we have 


p = -In (1 - R^jj) . 


Finally, t n is determined from the inequality 


T i T o T ___ I T, T, T 

-z— + + • . . + -Z \ p S -zr~ + -rr- + • • . + 

T 1 T 2 T m-1 T 1 T 2 T m 


and the equation 


C D = T 1 + T 2 + + T m-1 + T m [P - (T7 + tJ 


+ . . . + 


Sari 

T ra -1 


)] 


Once tp is determined, the daughter decay partial weight can be 
calculated in a more straightforward manner by simply summing up the partial 
exponential weights for the remaining intervals from T m up to T n . 


W, = exp [ 




t 2 

+ — + 
t 2 


T m (daughter) 


m 

T m 


- P) 


x m+l , , '•n l 

T m+ i (daughter) x n (daughter) ■> 


The trajectory calculation proceeds generally as before (Combi and Delsemme 
1980a, Combi 1980) except that the radiation pressure acceleration is also 
time dependent (through the heliocentric distance and velocity) in general and 
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must also be calculated in finite time steps. This calculation proceeds as 
the daughter weight terms are summed. 

A few preliminary model runs have been made with the new time dependent 
model. As one might expect, if one looks at densities or column densities 
reasonably close to the nucleus (£10^ km), there is little difference from a 
simple steady-state model. This first-order equivalence is due to the fact 
that as the comet changes heliocentric distance, photochemical lifetimes and 
parent molecule production rates change in opposite directions, tending to 
cancel each other. This cancellation also depends upon the production rate 
variation assumed in the model. If a simple inverse square law is assumed, as 
is true for small heliocentric distances (r H <0.7 AU) in the case of water 
vaporization or for the case of vaporization of something more volatile such 
as CO or CO 2 , the cancellation in the absence of radiation pressure is very 
strong. However, for a true water vaporization curve, the time dependence in 
the model becomes much more important. 

Figure 7 shows calculations for the CN coma as produced by HCN photodis- 
sociation with the new time dependent MCPTM at a heliocentric distance of 1 AU 
at times both pre- and post-perihelion. The orbit of Comet West was assumed 
as before; also, the vaporization curve for a dusty ^0 nucleus (Weissman and 
Kieffer 1984) was adopted. These model results at 1 AU imply an overall 
asymmetry of 10 percent in the range of 10^ to 10“* km usually covered by 
photometric apertures. The asymmetry increases to 20 percent for similar 
model runs at 2.5 AU. 

Perhaps the importance of time dependence is more in the interpretation 
of the limited spatial distribution data which are used to characterize the 
Haser scale lengths photometric observers typically use to reduce their data. 
Daughter scale length determinations are usually quite sensitive to the outer 
portion of a brightness profile, and this is just the region where time depen- 
dence (or really its usual neglect) is important. Figure 8 shows the sunward 
and antisunward brightness profiles for the same two model runs as above. As 
in the case for collisions, radiation pressure seems to enhance also the 
effect of time dependence on the CN spatial distribution, especially for the 
sunward profile. 
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Modeled CN(O-O) Intensity (Rayleighs) 


1 AU pre-perihelion 



1 AU post-perihelion 



Figure 7. CN(O-O) Intensity Maps for' Comet P/Halley at 1 AU 
Pre- and Post-Perihelion. These maps were 
computed with the new time dependent Monte Carlo 
particle-trajectory, model for the orbit of Comet 
P/Halley assuming a water vaporization curve and 
the same production rate for CN of lx 10^ s 
The asymmetry introduced by the time dependence is 
most noticeable in the sunward direction and 
decreasing as one approaches the antisunward 
direction. The same g-factor was assumed for both 
maps in order to illustrate the actual abundance 
asymmetry. 
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Modeled CN(O-O) Brightness Profiles 




log distance Ckm) 

Figure 8. Sunward and Antisunward Brightness Profiles for CN(O-O) in 
Comet P/Halley at 1 AU Pre- and „ost-Perihelion. As seen in 
Figure 7, the pre- to post-perihelion asymmetry is not present 
exactly antisunward but is quite noticeable in the sunward 
direction. See the caption for Figure 7 for modeling details. 
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